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Abstract 

Tai-Ping Liu |12j introduced the notion of “physical solution” of the isentropic Euler 
system when the gas is surrounded by vacuum. This notion can be interpreted by saying 
that the front is driven by a force resulting from a Holder singularity of the sound speed. 
We address the question of when this acceleration appears or when the front just move 
at constant velocity. 

We know from muz] that smooth isentropic flows with a non-accelerated front exist 
globally in time, for suitable initial data. In even space dimension, these solutions may 
persist for all t E M ; we say that they are eternal. We derive a sufficient condition 
in terms of the initial data, under which the boundary singularity must appear. As a 
consequence, we show that, in contrast to the even-dimensional case, eternal flows with a 
non-accelerated front don’t exist in odd space dimension. 

In one space dimension, we give a refined definition of physical solutions. We show 
that for a shock-free flow, their asymptotics as both ends t —>■ ±oo are intimately related 
to each other. 


1 The isentropic Euler equations 

We consider a compressible isentropic gas of density p, velocity u, pressnre p and specific 
internal energy e. The flow is governed by the Euler system 

(1) dtp + dw{pu) = 0, 

(2) dt{pu)+ Bw{pu®u)+ Vp{p) = 0. 

Classical solutions of ([TIE]) satisfy in addition the conservation of energy 

dt{^p\u\‘^ + pe) + div((^p|n|^ + pe + p)u) = 0. 
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When the solution develops singularities, say shock waves, the latter equality is replaced by the 
inequality 


(3) 



which holds true in the distributional sense ; Q plays the role of an admissibility condition. 
For the sake of clarity, we assume a polytropic equation of state 


p{p) = Ap^, A>0, 

where 7 > 1 is the adiabatic constant. Then the internal energy is given by 



(4) 


The sound speed is c = = \/'jA p^ where k := . We recall that a mono-atomic gas 

corresponds to the choice 



which means that the molecules have only d degrees of freedom, all of them associated with 
translations in space. In the physically relevant case d = 3, 73 = | is the parameter associated 
with mono-atomic gases like Argon. Di-atomic gases, like H 2 , O 2 or the air, obey to a pressure 


law with 7 = I < 73 . If d = 2, then 72 = 2 corresponds to the shallow water equations. 


We are interested in this paper in flows for which the total mass and energy are hnite, and 
the gas occupies a bounded region fl(t) of the ambient space at time t. The gas is surrounded 
by vacuum. We assume a priori that the front r(t) = dQ{t) is a smooth or piecewise smooth 
hypersurface in M'^. We denote 


n = {(x, t) I a: G fl(t) and t G M}, T = {(x, t) | a; G T{t) and t G M}. 


For weak admissible solutions of the Euler system, the total mass is conserved and the total 
energy is a non-increasing function of time : 



The latter inequality is an equality if the solution is classical. 

Of particular importance is the nature of the boundary condition along F(t). On the one 
hand, because there should not be any transfer of mass accross the boundary, F moves at the 
velocity u. More precisely, its normal velocity equals u ■ u, where u is the unit outer normal. 
The boundary can be viewed as a collection of particles moving at the fluid velocity u ; their 
paths are integrals curves of the ODE 


dX /X X 

— = u{X{t),t)- 


(5) 
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On the other hand, the Rankine-Hugoniot relations reduce on F to [p(p)] = 0, where p vanishes 
on the vacuum side ; therefore p must vanish on the interior side, meaning that p is continuous 
across T{t) (see also |II]). We mention in passing that this property might be violated at isolated 
times because an isolated (in space and time) discontinuity is not a shock front. An explicit 
example of that possibility was given by Greenspan & Butler [5], see Section[2]( “sub/supersonic” 
example). More generally, we allow discontinuities that occur on a codimension-1 subset of F, 
because they are not shock waves. 

Once we know that p = 0 along the boundary, there remains to understand its regularity, 
or singularity, in terms of the distance to the F(f). Here one can think of two approaches, 
complementary to each other. Both are based on the quasi-linear symmetric form of the system, 

(6) c)iC-I- (u ■ V)c -I- Kcdivu = 0, 

(7) dtu + {u ■ V)u + KcVc = 0, 

where c := - is a renormalized version of the sound speed. We point out that system (16117)) 
removes the singularity at p = 0, inherent to the Euler system, whose hyperbolicity degenerates 
at vacuum. These approaches distinguish two regimes, depending on whether the quantity cVc 
vanishes or not. If it does, then the particles at the front move freely ; we say that the front 
is not accelerated, or that the flow is smooth up to vacuum. If instead cVc ^ 0, then the front 
experiences a normal acceleration. This is what T.-P. Liu called a physical singularity at the 
boundary. Of course it requires that c be ^-Holder at the boundary. Beware that a given flow 
can be smooth up to vacuum in some region, for instance on some time interval, and display 
a physical singularity elsewhere. The transition between both regimes is still not completely 
understood. 

Outline of the paper. Section 12.21 deals with flows for which the sound speed is smooth 
enough at the vacuum, so that the front is not accelerated. Such flows arise for convenient 
initial data, but might exist only for a finite time. When the boundary regularity is lost, it 
is expected that a Holder-type singularity of c develops at the boundary. These singularities, 
called physical by T.-P. Liu, are investigated in Section [3l 

Our main results are Theorem 12.51 (non-existence in odd space dimension of eternal flows 
that are smooth up to vacuum) and Theorem 13.11 on the relation between the asymptotics as 
t ±cx3 in one space-dimension. 

2 Flows that are smooth up to vacuum 

The hrst approach considers a solution {c,u) of flMT)) . smooth in the entire space with a 
compactly supported component c. The following result is a classical application of the theory 
of symmetric hyperbolic systems, which can be found in ^ . 
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Theorem 2.1 Let the initial data (cq, 'Uq) belong to the Sobolev space for some param¬ 

eter s > 1 + I . Then there exists an open time interval I containing t = 0, and a unigue 
local-in-time solution in the class C{T, H^) r\C^{T, 

Thanks to Sobolev embedding, this solntion is classical : (c, u) G C^(/ x M'^). If cq is compactly 
snpported, then c(-,t) is too, and c vanishes at the bonndary of its snpport. 

The following argnment is adapted from Lin & Yang [13], where it was designed in the 
context of damped flows. For a solntion given by Theorem 12.11 the eqnation ([7]) rednces to 

(8) dtu{u ■'V)u = 0 along T. 

Since T is transported by the flow, this exactly means that the trajectories dehned by ([5]) and 
originating on the bonndary (that is, Y(0) G r(0)) remain on T and have a constant velocity 

u{X{t), t) = cst. 


We snmmarize this argnment into 

Proposition 2.1 Let {c,u) be a flow given by Theorem, \2.1\ with fl(t) bounded. Then the front 
r(t) is obtained from r(0) by a (uniform in time) transport: 

r(() = Mm) 

where 

fliflx) = X + tuo{x). 

We point ont that if is not the flow map of the gas, althongh it coincides with it on the vacnnm 
bonndary. 

Becanse they are classical ones, the solntions provided by Theorem 12.11 yield perfectly ad¬ 
missible solntions of the Enler eqnations ffT]|2|) . At this stage, it seems that we don’t need an 
extra condition at the bonndary. Their only flaw is that they are dehned only on a hnite time 
interval, bnt this is something we are accnstomate with in the theory of hyperbolic conserva¬ 
tion laws. We anticipate that if / = (—r*,T*) is bonnded for a maximal solntion, then some 
kind of singnlarity will appear as f —)■ T*. We have in mind the formation of shock waves as 
usnal, bnt becanse the front T = dfl is a place where the hyperbolicity of the Enler system 
degenerates, the singnlarity might appear at r(T*) and then propagate along T. We shall see 
later on snfhcient conditions that lead to snch bonndary singnlarities in hnite time. 

2.1 Global smooth solutions 

It may happen that T* be inhnite, in which case the smooth solntion of the forward Canchy 
problem is global-in-time. Snch solntions where hrst constrncted in ira when uo is close to a 
linear held x i—)■ Ax, and the spectrnm (the set of eigenvalnes) of the matrix A is contained in 
C \ (—oo, 0]. This resnlt was soon generalized to more realistic initial data : 
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Theorem 2.2 (M. Grassin [7].) In Theorem \2.1[ assume in addition that the spectrum of 
Vuo{x) is contained in some fixed compact subset of C \ (—cxo,0]. There exists an e > 0, 
depending upon uq, such that if 

||co||_ffo < e, 

then the solution exists for all positive time : T* = +cx). 

The idea behind the proof of Theorem 12.21 is that the data is close to (0, uq). For the latter 
data, the solution is (c = 0,n), with u governed by the vectorial Burgers equation 

(9) dtu + {u ■ V)u = 0. 

By assumption, the flow fifix) = x + f'Uo(x) is one-to-one at every positive time and therefore 
the solution of ([9]) is smooth for t > 0. For a non-zero Cq, there is a competition between the 
dispersion induced by -0* (the distance between two particles tends to increase linearly in time), 
and the nonlinear effect due to the pressure. Because of the dispersion, the density decays at 
an algebraic rate, and therefore the role of the pressure force gets smaller and smaller. If it 
was weak enough at initial time, we may expect that it will never be strong enough to lead 
to shock formation. An other important point is that the characteristic cones, which move at 
velocity u + with |(^| = c(p), will never overtake the boundary F(f), because the slope of c at 
the boundary is dominated by the diverging gradient of u. 

Eternal smooth flows 

Theorem 12.21 has the beautiful consequence that smooth eternal compactly supported flows 
exist in even space dimension : 

Corollary 2.1 (d even.) Suppose the space dimension d is even. Then in Theorem \2.1[ one 
can choose the initial data so that the solution is eternal .■ / = M. 

Proof 

Just choose Uq so that the spectrum of Vmo(x) is contained in a hxed compact subset K of 
C \ R. Then choose cq small enough in Because K does not meet (—cxo,0], we have 

T* = -l-cxD. Because K does not meet [0, J-oo), we have T* = —oo. 


Corollary 12.11 raises the question of whether there exist eternal solutions with compact 
support, smooth up to vacuum, in odd space dimension. The construction in the proof above is 
not possible, because odd-size matrices do have at least one real eigenvalue. It is remarkable that 
one meets an obstruction, see Theorem 12.51 below. The non-existence in one space dimension 
is a consequence of a calculus due to P. Lax ra : the Euler system, whose wave velocities are 
X± := u ± c, can be diagonalized in Riemann coordinates 

r± = u ± c, 
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namely (I6]I7|) is equivalent to 


{dt + \±d^)r± = 0. 


Each of the transport equations above tells us that r± keeps a constant value along every char¬ 
acteristic curve t X{t) dehned by ([5]). Lax showed that there exists a positive function N{p) 
such that the expressions y± := N{p)dxr± satisfy Ricatti equations along the characteristics 


r)\ 

(dt + X±dx)y± + N = 0. 

OTj- 


( 10 ) 


Because of genuine nonlinearity 


the only eternal solution of flTUD is = 0. Therefore the global existence implies that r±, 
or equivalently p and u, are constant. With the assumption of hnite mass, we conclude that 
p = 0 : there is no fluid at all. In the next section, we shall use the fact that if a solution y of 
oa exists on a time interval (T, -|-cx)) (respectively on (—oo,T)), then p > 0, that is dxr± > 0 
(resp. dxr± < 0 ). 

Notice that if 7 = 71 = 3, the system above decouples as two independent Burgers equations, 
since then r± = A± : 


( 11 ) 


{dt + X±dx)X± — 0 . 


In this case, one can take = 1. 

We summarize this analysis into 

Proposition 2.2 {d=l.) In one space dimension, there does not exist a non-trivial eternal 
solution (c, u) of m. of class C^, with c of compact support. 

2.2 Dispersion vs smooth boundary 

The following statement tells us that if the initial velocity held does not drive fast enough the 
boundary, then a singularity must happen in hnite time. 

Theorem 2.3 We assume 7 < 1 -|- . 

Let (co,Mo) £ he an initial data, with s > 1 + f ■ Assume that cq is compactly 

supported in ff(0), a domain with smooth boundary. Assume that 


( 12 ) 


.2(0) 

Then the maximal time T* > 0 of existence provided by Theorem \2.1\ is finite: 



det Vuo{x) dx < 0 . 


T* < -Fcx). 
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Remark that the integral J in flT^ depends only upon the restriction of uq to the boundary 
(9f2(0). The identity 

det Vtto(a;) = div(MQVMo x Vtig) 

yields the alternate formula 

/ det Vmo(3;) = / uldet{u,'\/ul,'\/ul) ds{x), 

Jn{o) Jdn{o) 

which shows that J depends on Uq only through its restriction to the boundary. By symmetry, 
the same is true for every component Mq. 

To prove Theorem 12.31 we begin with a property which must be rather classical. A simpler 
version, can be found in Milnor’s paper [16] . Notice that this lemma is valid for every parameter 
7 > 1. 

Lemma 2.1 As a function of time, the volume Q{t) is a polynomial of degree < d. 


Proof 

A function G{t) is polynomial if and only if it is locally a polynomial. It is therefore enough 
to prove that G is a polynomial of degree < d over some non-trivial interval (—e, e). 

Let us consider the map = x + tuo{x). Whenever |t|||uo||i,oo is less than 1, %jjt is a. 

diffeomorphism. Therefore ?/^t(f2(0)) is the bounded open domain whose boundary is (512(0)). 
Because 'ipt coincides with the flow of the fluid over 512(0), we hnd that '0i(12(O)) is the bounded 
open domain whose boundary is 512(f), and we conclude 

fi(f) = V^i(fi(o)). 


We infer the formula 


|12(f)| = / detCViftix)) dx = / det{Id + tVuo) dx, 

Jq{o) Jq{o) 


where the integrand is a polynomial in t of degree at most d. 


Let us point out that the integral J is precisely the coefficient of the monomial in this 
polynomial. We now prove our theorem. 

Proof 

Let Q{t) be the polynomial 1 1 —)■ |12(f)|, Q{t) = + l.o.t. Our assumption amounts to 

(13) |12(f)| <G((l + |f|)''-'), Vf>0. 
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Notice that if J < 0, we actually have Q{t) < 0 for large enough time, whence an immediate 
contradiction. However, it is important for the sequel to treat also the case J = 0. 

Let us recall the following inequality, which can be found in mm- It is a direct consequence 
of the Euler system ([T1I2]) and of the energy inequality ([3]) : 


(14) 


d 

dt 


■^p\tu — xy + rpe ] dx <2{1 —dK)t / pedx. 


Remark that the quantity in the left integral is a linear combination of the density, the mo¬ 
mentum and the mechanical energy with polynomial coefficients : 


1 , 


x\^ -I- t^pe = 


+ X • {pu) + f + pe 


When 7 G (1, 7 ^], the factor 1 — dK is non-negative. With the Gronwall inequality, we deduce 
the estimate 


(15) 



p'^dx = cst 


pedx = O ((1 -I- t) . 


Now, we can estimate the total mass with the help of Holder Inequality: 


M= pdx < 
Jn{t) 


\ 1/7 
p'^dxj 


= O (^(1 + j = O (^(1 + . 

Letting t —>■ -|-oo, we deduce that M = 0. In other words, there is no fluid at all. 


If instead 7 > 7 d, then 1 — dK < 0 and we have 


r' + CxD 


tdt pedx < 


dn — 1 


\x\^Pq{x) dx, 


wich gives 

j‘+oo 

/ tdt j p^dx < 00 . 
Jo Jr'^ 

Using the same Holder inequality as above, we infer 


“+00 




f(l + < CX), 


which implies either M = 0 (no fluid at all) or 


1 + (d- 1)(1 - 7 ) < - 1 , 





Q.E.D. 


that is 7 > 1 + . 


Remark that Theorem 12.31 doesn’t fully use that (c, u) is a classical solution. We only need 
on the one hand that (p, pu) be a weak solution of the Euler system, satisfying the “entropy” 
inequality (j3]) (this is admissibility in the sense of Lax), and on the other hand that (c, m) is of 
class in a neighborhood of the front T. There is a statement analogous to Theorem 12.31 in 
this context : 

Theorem 12.31 We assume 7 < 1 + . 

Let (p, pu) be a locally bounded admissible (in the sense that solution of the Euler system 
nViE) over (0, + 00 ) X M'^. Assume that p is supported in hi, a domain bounded in space at each 
time, whose boundary is smooth. Assume finally that (c, u) is of class W in a neighborhood of 
dfl. Then necessarily 

(16) / det Vmo( 2 ^) dx > 0 . 

Jn(Q) 

This shows that the development of shock waves is not sufficient to resolve the lack of 
smooth solutions. Many flows actually exhibit also some boundary singularity after some time. 
This is where T.-P. Liu’s notion of physical singularity comes into play. 


The mono-atomic case 

If 7 = 7 rf, Theorem 12.31 can be improved in quantitative way. Suppose that J > 0. From 

T 


d 

dt 


p\tu — xf + rpe dx <0, 


we infer 


A / p'^dx = / {'y — l)pedx < 


'n{t) 


in{t) 


7-1 


/, /:= 


'n(o) 


-po\x\‘^dx. 


Applying again the Hblder Inequality, this yields 

AM^ 


7-1 


f < |f2(t)|2/'’*/. 


With |f2(f)| = Jt‘^ + l.o.t., we deduce that for a solution to be smooth at the boundary for all 
f > 0 , we must have 

AM'^ 9/^ 


7-1 


We therefore have 
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Theorem 2.4 We assume 7 = 1 + . 

Let (co,Mo) e H^{W^) be an initial data, with s > 1 + |. Assume that the domain 

fi(0) = {co > 0} 

is bounded, with smooth boundary. Assume that 

(17) J := / det Vtto(a^) < 

Jn{o) 

Then the maximal time T* > 0 of existence provided by Theorem A2. 1\ is finite: 

T* < +CX). 


2AM^ 


The odd-dimensional case 

As anounced above, we prove that an obstruction occurs in the much more general context of 
an odd space dimension, a case which includes the realistic 3-dimensional space. The following 
result answers a question raised in [18] about eternal solutions. 

Theorem 2.5 {d odd.) We assume that the space dimension d is odd, and that 7 < 1 -|- . 

Then there does not exist a non-trivial eternal solution {c,u) G of (0[^ (hence 

smooth up to vacuum), with cq supported in a compact domain with smooth boundary. 


Proof 

Suppose (c,u) is such a solution, and denote Q{t) = |fl(t)|. From Lemma EH] Q is a 
polynomial of degree < d. Because Q takes positive values for every f G M, we hnd that deg Q 
must be an even number. Because d is odd, we deduce that actually deg Q < d — 1, whence 

(18) |S1 (()|=0((1 + |(|)''-'). 

The rest of the proof is exactly the same as in the proof of Theorem 12.31 

Q.E.D. 

Again, Theorem 12.51 doesn’t fully use that (c, u) is a classical solution, but only that it 
is smooth up to the vacuum. It tells us more about the onset of physical singularity at the 
boundary, which we consider in the next Section. 
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3 Physical singularity 


It was first recognized by T.-P. Liu jT2] that many flows that are classical solutions in the 
interior of Q{t) must experience a singularity at the boundary. His hrst motivation was the fact 
that for the Euler system with damping, the front travels so slowly that Q{t) remains uniformly 
bounded. For gaseous stars, the gravity has a similar effect [15] . This lack of dispersion comes 
in conflict with the decay of the integral of and the conservation of the total mass. We have 
shown in Theorem l2.3f that such an obstruction is present even without damping. Since the gas 
must flow anyway, the way to resolve the obstruction is too admit that c is only Holder-regular, 
of exponent |, at the boundary. In other words, is Lipschitz and the normal derivative does 
not necessarily vanish 


(19) 


dc^ 


9 ■= 


-y — 1 du 


> 0 . 


That g is non-negative follows from the fact that is positive in f2(t) and vanishes on the 
boundary. Assuming that u is smooth accross the boundary, we obtain the identity 


^u{X{t),t) = gu 

along a boundary path. Therefore g can be viewed as an acceleration of the front expanding 
in vaccum. The possibility that g be non-zero resolves of course the obstruction raised by 
Theorem 12.31 . It allows the volume of the domain to grow fast enough, at least as whatever 
the initial velocity. Notice that we don’t need that g be positive for every time. For instance, 
an appropriate initial data in even dimension yields a global smooth solution (Theorem 12. Ih . 
for which we just have g = Q. 

In one space-dimension, the local existence of a flow with physical singularity at vacuum 
was proved in [IT] and [2] by two different methods. These works deal with the regime where 
g is strictly positive ; it seems to be an open problem to have an existence result covering both 
regimes g = ^ and 5 ^ > 0, and in particular the transition from one regime to the other. The 
local existence in several space dimensions is proved in [HIE]- 

Global existence of a weak entropy solution in one space dimension, in presence of vacuum, 
has been proved by means of Compensated Compactness, see mu- However, the method of 
proof says nothing about the behaviour of the solution at the front with vacuum. At least, 
this approach has the merit to deal with both regimes, the accelarated and the non-accelerated 
ones. One drawback is that the flow is obtained as the limit of approximated solutions, whose 
support grows faster than that expected for the genuine flow. 


3.1 One-dimensional physical singularity 

We wish to give an accurate definition of what is an admissible flow, from the point of view of 
physical singularity, at least in the one-dimensional case. Let Q{t) = (a(f), b{t)) be the domain 
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occupied by the fluid (where p(-,t) is positive)0. We have b'{t) = u{b{t),t) and therefore 
b" = g > 0 ] hence b is convex. Likewise a is a concave function. Notice that the discontinuity 
points of a' or b' form sets that are either finite or denumerable. 

Let us begin with the mono-atomic case, where 7 = 8 . We have already seen that the 
corresponding Euler system splits into two coupled Burgers equation. Therefore, as long as 
the flow is smooth, the characteristics are straight lines. Even more, the tangent lines to the 
boundary are characteristics. The tangent at a point P actually splits into two halves, one being 
a 1-characteristic and the other a 2-characteristic. This is clear when considering a point Q 
close to the boundary, say to the right component T^ : because b is convex, there are exactly two 
tangents to T^ passing through Q, which are the 1- and the 2-characteristics respectively. We 
observe that a 2 -characteristic cannot emanate from T^, and a 1 -characteristic cannot terminate 
in T^. 

When 7 > 1 is a general parameter, our definition of an admissible flow is that the same 
pattern be true. 

Definition 3.1 {d = 1.) Let 7 6e > 1. An admissible flow surrounded by vacuum is a measur¬ 
able hounded field (p > 0 ,m), which satisfies the following requirements: 

1. The domain Vt is hounded at left and right by two curves t a{t),b{t), with —a and b 
convex, 

2. The flow is a distributional solution of /[T[E) in 12 , 

3. It satisfies the energy inequality 

4- For almost every boundary point P G dVL, 

lim pix, t) = 0, 

(x,t)^P 

5. If a 2-characteristic fi 2 reaches T,, at some point P, then P is the terminal point of 
h ?■ if ^2 reaches at some point Q, then Q is its initial point. Symmetrically, if a 
1-characteristic fii reaches T^ at some point P, then P is the initial point of fii ; if fii 
reaches at some point Q, then Q is its terminal point. 

The reader will have noticed that the definition above is somewhat sloppy : without some 
regularity, characteristic curves may not be well-defined. We shall therefore use it only in 
situations where (p, u) are smooth enough. This will be the case in Theorem 13.11 Notice 
that our definition is consistent with the rarefaction waves displayed in the “sub/supersonic” 
example below ; this is a case where infinitely many characteristics emanate from the same 
boundary point. 

^We do not consider the case where the number of connected components varies with time. 
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Examples 

We give below two explicit examples of eternal flows with physical singularity at vacuum. They 
are built in one-space dimension for a mono-atomic gas (7 = 7 i = 3), for which the Euler system 
decouples as a pair of Burgers equations away from shock waves. Without loss of generality, 
we set so that c = p. The wave velocities \± = u ± p satisfy 

( 20 ) dtX + Xd^X = 0 . 

In a domain where A is Lipschitz, the characteristic curves, on which A is constant, are lines of 
slope A. 

We point out that in both examples, uq vanishes identically. Because our solutions are 
smooth in hi, this implies the reversibility : 

p{x,—t) = p{x,t), u{x,—t) = —u{x,t). 

An accelerated case. In the following example, the front between the gas and vacuum is the 
hyperbola dehned by the equation 


x"^ = 1 + t^. 

In the gas, the density and velocity of the flow are given by 

tx 

We leave the reader verifying that each of the functions 

tx ± vT+l^"--^ 


X±ix,t) = 

satishes the Burgers equation. Obviously, 


l + f2 


c2 = 


1 + r -x^ 


(l+t2)2 

vanishes at the boundary, where it is Lipschitz, and the acceleration 

1 


g{y/l + t 2 ,t) = 


(1 + f2)3/2 


is positive for all time. Mind however that this expression is integrable in time. The 
domain Q{t) behaves asymptotically as the interval {—\t\, |t|), and the flow is asymptotic 
to 




X 


pR{x,t) = -^1- — , UR{x,t) = -. 
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The physical singularity is damped as |t| —)■ oo ; the acceleration g decays like \t\ 

Let us point out that this example does not contradict Proposition 12.21 because even 
after extending c by zero outside of hi, the regularity of the solution (c, u) is too low to 
implement Lax’ calculation. As a matter of fact, the Cauchy-Lipchitz theorem does not 
apply to the characteristic flow because the wave velocities are not Lipschitz : P is the 
envelop of the characteristic lines. 

A discontinuous example. This one is mimics a situation described in | 8 ]. Let us choose 
instead the initial data 


po = 1, Uo = 0 in f2(0) = (-1,1). 

The solution (c, u) of fEH7|) is Lipschitz, except at the points {x,t) = (±1,0) where a 
discontinuity happens (obvious in the data above). At time t > 0, the gas occupies the 
domain Q{t) = {—t — 1, f ± 1) and the solution is given by 

f if —t — 1 < a; < t — 1, f if —t — 1 < a: < —f ± 1, 

if f — i^x<t + l. ^ ^ I ^ if —f±l<a:<t±l. 

The flow is a rarefaction wave in the domain defined by |t — 1| <a;<f±l, as well as in 
that defined by —f — 1 < x < —|f — 1|. The front is accelerated only at time t = 0, but 
this acceleration is an impulse : the velocity of the front flips from —1 to ±1. 


3.2 The domain in one space dimension 

Let us consider an eternal flow surrounded by vacuum in the sense of definition 13.11 We suppose 
in addition that p,u are smooth in hi. We are interested in the growth of |fl(t)| as t ^ ± 00 . 
Because of the decay of the integral of p^, we already know that |fl(t)| is bounded below by 
for some constant Ci > 0 . 

For a one-dimensional gas, the domain Q{t) = {a{t),b{t)) is an interval where t ^ —a,b 
are convex functions. The derivative b' has limits q- < q+ as t ^ ± 00 . Likewise a' has limits 
p+ < P-. Because of the lower bound of the volume, we know that 

P+<(1+, q-<p-. 


We actually have 

Theorem 3.1 {d = 1.) Let a flow be smooth (shock-free) in its domain Vl{t), where 12(0) = 
(ao,feo) is a bounded interval. We assume that the flow is admissible in the sense of Definition 

Then p+ = q_ and p_ = . 
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Proof 

We consider 2-characteristics, which are integral curves of the ODE 

dX w N 

-— = X{X,t), X := u + c. 

(JjL 

Because the flow is smooth in D, we know that r := w + f is constant along each 2-characteristic. 

We choose the origin of the time arrow in such a way that a' and h' be continuous at t = 0. 
Denote D"*" = D 0 {t > 0}. By Cauchy-Lipschitz, D’*' is foliated by the 2-characteristics. On 
such a curve /3, the time goes from > 0 to At t = Ti^ (3 reaches the boundary 

at some point m. By admissibility, we actually have 


m e r+ := r+ U (ao, 6o), r+ := w n {t > 0}. 


We order the set of 2-characteristics in D’*' from left to right, and denote it (—cx),/9o), where 
/So = {(^ 0 , 0 )} and —oo is the “point at inhnity” along B^. 

Consider two such characteristics (5 -< f3'. If TP^{(3) is hnite, the terminal point /3(T'f*"') is 
on r^. Then /3' is contained in the bounded set delimited by (3 at left and T^ at right, and we 
infer TP"'{/3') < TP'^{/3) < oo. We deduce that f3 i-)- PP^ is non-increasing, and there exists 
(3* G [—cxD,/?o] such that 

^ (T/-(/3) = +oo), 

(/S*^/3) ^ (T^“(/5) <+oo). 


In the limit case where pP'^ is hnite for every (3 (unlikely), then j3* = —oo. If instead T-f*" = -|-cxd 
for every (3 (unlikely), then j3* = h^. 

When 13 -< 13*, the Ricatti equation along f3, plus the fact that 


dr 


7 + 1 
4 


> 0 , 


imply that dxr > 0. This shows that (3 e-)■ r\p is non-decreasing up to (3*. On the other hand, 
if j3* + (3, we have seen that (3 i—)■ pP^((3) is non-increasing. Since r equals u = b' over the 
boundary, and b' is non-decreasing, we deduce that /3 i—)■ r |/3 is non-increasing beyond /3*. If 
/3* -< bo, we have 


= sup u < supr < supr = r|/ 3 » = sup r |/3 < sup u = q+. 
r+ r+ 4*^4 r+ 

If instead (3* = bo, then 

qj^ = sup-u < supr < supr = lim r |/3 = lim r{x,0) = r{bo,0) = b'{0) < g+, 

P+ P+ ^ ^ ^^0 

where the last equality stands because b' is continuous at t = 0. In all cases, we infer 

sup r = r|/ 3 * = g+. 
o+ 
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Considering symmetrically the foliation of fl {t < 0} by the 2-characteristics, we 

hnd a /3* G [ctO) + 00 ] snch that 


^ (r-(/ 5 ) >-00), 

(r"(/ 3 ) =-00). 

We have that (3 1 —)■ r|^ is non-decreasing over (oq, (3*), non-increasing over (/?*, -|-oo), and satishes 

snpr = r\/s^ = p_. 
n- 

We point ont that those 2-characteristics passing throngh a point m G (ao,foo) belong to 
both sets, associated with either or fl_. 

Snppose now that q+. Withont loss of generality, we have p_ < q+. Then r|^» > r|^ for 
every 2-characteristic (3 in Thns (3* does not emanate from a point m G (oq, bo) ; we mnst 
have (3* -< ao. Thns (3* emanates from some point of T^ and we have r|^* = r{a{t),t). 

This gives q+ = a'{t) < p_, a contradiction. 

This ends the proof of the Theorem. 


Open problems 

One important issne is of conrse to give a mnlti-dimensional version of Dehnition l3.ll The difh- 
cnlty here is that we cannot distingnish between two characteristic families : the characteristic 
cnrves passing throngh a point P are tangent to a cone, which is a connected set. This cone 
shrinks to a line at bonndary points. 

Becanse of the dispersion estimate and the conservation of mass, we know that |t|'^ = 
0(|D(t)|) as |t| —)■ -fcxD. We snspect that the converse ineqnality |D(t)| = is also trne. 

Theorem 13.11 provides a positive answer in the special case where d = 1 and the how is shock- 
free. It is also trne for hows that are smooth np to the bonndary, see Lemma 12.11 

Actnally, in both cases, the domain VL{t) behaves asymptotically as tO for the same O as 
t —)■ ± 00 . This latter property is donbtfnl for hows that display shocks, bnt is it trne for 
shock-free hows with physical singnlarity when d > 2 ? 
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